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Abstract 



The Wnt signalling pathway transducing the stabilization of /3-catenin is 
essential for metazoan embryo development and is misregulated in many dis- 
eases such as cancers. In recent years models have been proposed for the 
Wnt signalling pathway during the segmentation process in developing em- 
bryos. Many of these include negative feedback loops build around Axin2 
regulation. In this article we propose a new negative feedback model for the 
Wnt pathway with Dickkopfl (Dkkl) at its core. Dkkl is a negative regula- 
tor of Wnt signalling. In chicken and mouse embryos there is a gradient of 
Wnt in the presomitic mesoderm (PSM) decreasing from the posterior to the 
anterior end. The formation of somites and the oscillations of Wnt target 
genes are controlled by this gradient. Here we incorporate a Wnt gradient 
and show that synchronization of neighbouring cells in the PSM is important 
in accordance with experimental observations. 

Keywords: Clock and wavefront model / Dkkl / Modelling / Segmenta- 
tion / Wnt pathway 

Introduction 

A segmented body plan is a fundamental characteristic feature of vertebrates. 
The process of segmentation is carried out by a combination of changes in 
gene expression and relative anterior-posterior cell position in the presomitic 
mesoderm (PSM) (?). In the anterior end of the embryo the somites are 
segmented at equally separated time points with species dependent periods. 
In mice the period is around 120 min and in frogs it is around 90 min. 

In (?) Cooke and Zeeman proposed the clock and wavefront model to 
describe the segmentation process. The idea is that locally coupled oscillators 
are controlled by a morphogen gradient in the PSM. The oscillators are the 
clocks providing temporal information, e.g., cycle state, and the morphogen 
gradient is the wavefront providing spatial information about axial position. 
Until now three major pathways controlling the segmentation process have 
been found: the Notch, Wnt and FGF pathways. They all have target genes, 
which oscillates and, interestingly, Wnt target genes oscillate out of phase 
with Notch and FGF target genes (?). These three pathways could be the 
clocks. There is a decreasing gradient of wnt3a and fibroblast growth factor 
8 (fgf8) starting from the tail bud through the PSM (??). Fgf8 seems to be 
regulated by Wnt (?), thus Wnt could form the main wavefront. The actual 
setting of the somites happens at the determination front, where the gradient 
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drops below a certain threshold. Cells past the determination front become 
permissive to form somites depending on their phase of oscillation (??). 

In 2003 it was discovered by Aulehla et al. that Axin2 oscillates during 
the segmentation process in developing mouse embryos. Since their discovery 
several models for the Wnt oscillator have been proposed (???) with Axin2 
as the main variable. We propose a new model for the Wnt pathway with 
Dickkopfl (Dkkl) as the core of a negative feedback loop. The reason for 
changing from Axin2 to Dkkl is that Axin2 deleted mice (Axin2~/~) have 
no dramatic phenotype (?). Dkkl has an oscillatory behavior during the 
segmentation process in mouse embryos (?) and lack of Dkkl in mouse 
embryos results in anterior head truncations and neonatal death (?). In mice, 
for which we want to model the segmentation process, several Dickkorpf genes 
are known. Dkkl is the most widely studied of these. 

Results and Discussion 
Modeling the Wnt//3-catenin pathway 

During Wnt signaling /3-catenin interacts with the TCF/LEF-1 DNA-binding 
proteins to promote transcription of Wnt target genes (?). As for Axin2 
the transcription factor for Dkkl is /9-catenin (??). After transcription and 
translation Dkkl goes through the cellular membrane where it can bind to the 
extracellular domains of the low-density lipoprotein receptor-related protein 
5 and 6 (LRP5/6). When bound to LRP5/6, Dkkl acts as an inhibitor of 
Wnt signaling by blocking the association between Wnt, Frizzled (Fz) and 
LRP5/6 (?). Wnt acts as an inducer for the formation of this complex and 
Dkkl is a competitor to this induction (??). 

It has been proposed that the Wnt signal is transduced through the cell 
membrane by the binding of Dishevelled (Dsh) to the intracellular domain of 
the Fz receptor(?). Axin and Dsh can bind together via their DIX domains 
(?) and they co-localize at the membrane (?) during Wnt signaling. There- 
fore Dsh bound to Fz may recruit Axin bound to GSK3/3 to the LRP5/6 
receptor (?), where a phosphorylation of LRP5/6 is initiated. The LRP5/6 
receptor has a binding site for Axin and upon Wnt signalling GSK3/3 (bound 
to Axin) phosphorylates LRP5/6, which requires Axin (?). The phosphory- 
lated LRP5/6 receptor may be able to recruit and more efficiently bind the 
Axin-GSK3/3 complex to the membrane and the phosphorylation process is 
thereby amplified (?). 

At the cell membrane Axin is phosphorylated by GSK3/3 and then de- 
graded (??)• The degradation of Axin leads to a decrease in the formation 
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of the destruction complex comprised of /3-catenin, the kinases glycogen syn- 
thase kinase 3 (GSK3/3), casein kinase (CKIa), the scaffolding proteins Axin 
and adenomatous polyposis coli (APC). In the destruction complex /3-catenin 
gets phosphorylated and subsequently degraded. 

Interestingly enough GSK3/3 plays a dual role in controlling the Wnt sig- 
nal. When the Wnt signal is off then GSK3/3 phosphorylates /3-catenin in the 
destruction complex and when the signal is on then it phosphorylates Axin 
at the LRP5/6 receptor. Whether it is the same or distinct Axin-GSK3/3 
complexes that carry out the phosphorylation of /3-catenin and Axin is un- 
known (?). 

The dynamics of CKIa and APC are included in the parameters govern- 
ing the destruction complex, the TCF/LEF-1 dynamics are contained in the 
transcription of Dkkl, and the dynamics of Dsh are included in the formation 
of the complex consisting of Axin, GSK3/3 and LRP5/6 at the cell membrane. 
Figure [1] shows a simplified diagram of the proposed model and the associ- 
ated equations. C, [GA], G, B, L, D m , D, [LD], A and [LCM] are the 
concentrations of the destruction complex, GSK3/3-Axin complex, GSK3/3, 
/3-catenin, LRP5/6, Dkkl mRNA, Dkkl protein, Dkkl-LRP5/6 complex, 
Axin and LRP5/6-Axin-GSK3/3 complex. The formation and breaking of a 
complex X are denoted by Cfx and c&x, respectively. The transcription and 
translation rate of Dkkl are given by the parameters c tsc and c ts i, respectively. 
The Hill coefficient on B regarding the transcription of D m is associated to 
the amount of cooperativity between /3-catenin and the TCF/LEF-1 com- 
plex. For example, no cooperativity would result in a Hill coefficient of one. 

The concentration of GSK3/3 has been shown to be extremely stable (?) 
and consequently its total concentration, GSK3(3 tot , is assumed to be con- 
stant during the involved time scales. The same assumption goes for the total 
concentration of LRP5/6, since the half-life of LRP6 is around 4.7 hours (?). 
Therefore we have not included any source or sink for the concentrations of 
G and L. Only constitutive sources of /3-catenin and Axin, respectively, are 
included in the model, because free (unphosphorylated) /3-catenin and Axin 
is stable (?). 



Determination of the parameter values 

In ? a model for the Wnt pathway in Xenopus is proposed and from this 
article estimates for the parameters K c = > K[ga] — Sjrfj"> a ' GSK3/3 to (, 
Sa and Sb are taken. Dkkl binds to LRP5/6 with a high affinity, since the 
dissociation constant has been measured to be around -^[ld] = 0.4 — 0.5 nM 
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The other parameters are estimated to produce oscillations with a period 
of around 120 min and a very low concentration of total Axin as found by ?. 
The low concentration of total Axin is thought to act as a buffer to changes 
in the concentration of the other constituents, which may also take part in 
other signalling pathways (?). 

In the activation of transcription of Dkkl by /3-catenin we assume cooper- 
atively between /3-catenin and TCF/LEF. The Hill coefficient is set to three. 
A Hill coefficient of one and two returns no sustained oscillations, when the 
model is simulated. In other models of the Wnt pathway Hill coefficients of 
values 2 and 5 have been applied (??). 

Oscillations of the mRNA, protein and complex levels 

The parameter values produce oscillations of the involved constituents with 
a period of around 120 min as seen in Fig. [TJ In addition, the concentration 
of Axin is very low compared to the concentrations of the other variables, 
as discussed above. The different steps in the proposed Wnt model com- 
ply with the phase shifts in the oscillations of the different components. An 
increase in the [LCM.] concentration leads to a decrease in the Axin concentra- 
tion. This decrease causes a reduction in the concentration of the destruction 
complex and consequently an increase in the /3-catenin concentration. This 
increase will after a while cause an increase in the Dkkl concentration. The 
concentrations of GSK3/3 and [GA] are mirrors of each other, since a high 
concentration of will cause [GA] to leave less free GSK3/3 behind. 

? found no significant oscillations in the level of /3-catenin. For our 
choice of parameters /3-catenin show oscillatory behavior with an amplitude 
of approximately 5nM, which is not significantly low. A different set of 
parameters could possibly give a smaller amplitude of /3-catenin but the 
general results, presented later, is not significantly altered by this. 

A spatial Wnt gradient induced by time variation 

The described model could be the clock of the clock and wavefront model 
as mentioned above. The wavefront we introduce as a time dependent Wnt 
gradient through a time dependence of the parameter c/[lga]- The reason 
for making cj[lca] time dependent is clear if we introduce the variable [LW] 
describing the binding of Wnt (W) to the LRP5/6 receptor together with 
Frizzled. The temporal change of [LW] is described as 



= c f[LW] LW{t) - c b[LW] [LW] , 



(1) 



and if we assume steady state for the binding of Wnt to LRP5/6 then 



[LW] 



c fLW LW(t) _ LWjt) 

Cb[LW] K[LW] 



(2) 



Inserting this L estimate into the concentration of [LCM.] , it can be seen that 
the formation of [LCM.] will be time dependent, with 



To include the decreasing gradient of Wnt observed in mouse embryos we 
use a decreasing value of Cf[LGA]- If we assume that Wnt primarily decreases 
through a diffusive process, then the profile of cj[lga] will be Gaussian. In 
Fig. [3] the amplitude and period of Dkkl concentration is plotted as c^lqa] 
changes. The green line refers to the reference value of Cf[LGA]', reducing 
it results in smaller amplitudes. The period seems to be almost unaltered, 
which complies with the findings of constant periods of the oscillations in the 



The length of the PSM is approximately constant during the formation 
of the first 15-20 somites in mice embryos (?). The same is almost true for 
the size of the somites. From ? one can find that a cell budded off in the tail 
bud will be segmented within the order of ~ 1100 — 1300 mins. 

It has been measured that FGF (regulated by Wnt) exhibits a gradient 
in the PSM with a fold change of 2-5 (?). The fold change used for c^lqa] 
is 2. Assuming that Wnt diffuses through the PSM and setting the final 
value of Cf[LGA] equal to 2.5 enable us to calculate a Gaussian profile of 
Cf[LGA] representing the Wnt gradient. In Supplementary Fig. 1, Gaussian 
profiles of Cf[LGA] are plotted with different initial values. A decreasing value 
of Cf\iGA] i n the PSM will give rise to smaller amplitudes and slightly longer 
periods. 

Synchronization of neighbouring cells 

The elongation of the embryo is modelled by simulating cells budding off 
from the tail bud at regular time intervals, R. Elongation is considered in 
only one dimension. In Fig. |2]time series for cells with R = 10 min are shown. 
The final level of Dkkl at the determination front is oscillating (see Fig. [2]A.) 
with a period of 120 min. However, these oscillations are only possible if the 
onset of oscillations in the cells are synchronized, i.e. if the state for cell 
i at time t is denoted Si(t), then £1+1(0) = Si(R). This means that the 
initial state of a cell budded off from the tail bud equals that of its anterior 



Cf[LGA](t) 



Cf[LGA]W(t) 
K[LW] 



(3) 



PSM. 
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neighbour. If all the cells have the same initial state then the levels of Dkkl 
at the determination front will also be equal, see Fig. [2B. 

One could imagine that the synchronization of the clocks in each cell is 
not perfect. If the initial state of cell i + 1 is chosen randomly as Si+i(0) = 
Si(R + U( — ^, 7^)), then the period is almost unaltered (see Supplementary 
Fig. 2A). Here U(—x,x) denotes random numbers uniformly distributed 
between — x and x. If the initial state is randomly chosen within the range 
of Dkkl levels, then the oscillations of Dkkl at the determination front are 
disrupted and no periodicity is visual (see Supplementary Fig. 2B). Thus, 
the system seems to be robust to small changes in the synchronization. 

Diffusion of Wnt is thought to be the main reason for the Wnt gradient 
in the PSM. Thus a Gaussian profile of c^lga] is used. The profile of Cf^GA] 
would still be Gaussian if we also included a half-life of Wnt. If the Wnt 
gradient was only controlled by the half-life of Wnt then the levels of Dkkl 
still oscillate at the determination front with a period of 120 min, when the 
model is simulated as above with R = 10 min. 

The period of somite formation increases 

It is known that the period of somite formation increases as the embryo 
elongates (?). This led ? to propose that the level of Wnt production in 
the tail bud increases with time. To test their hypothesis we let the initial 
value of Cf[LGA\ increase with time keeping the level at the determination 
front fixed at 2.5 ~m^-=~ as above. Simulating an increase in the Wnt level of 
the tail bud results in longer periods of the Dkkl levels at the determination 
front, see Supplementary Fig. 3. Thus, we can confirm the hypothesis of 
?. Note that a steeper increase of the Wnt level in the tail bud results in a 
greater extension of the periods (compare Supplementary Figs. 3A and 3B). 

Outlook 

We propose a Wnt model with Dkkl as the core for the negative feedback 
loop, which exhibits sustained oscillations of Dkkl. By simulation of the 
elongating embryo we were able to test the importance of synchronization 
between neighbouring cells. In addition, we could also show that a small 
change in the synchronization did not significantly disrupt smooth oscilla- 
tions of Dkkl levels at the determination front. The clock and wavefront 
model was investigated through a decreasing Wnt gradient in the PSM and 
we could confirm the hypothesis that the Wnt level in the tail bud increases 
over time resulting in longer periods of the oscillations. Further aspects of 
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the growing embryo could be incorporated into the model, like a varying 
PSM size and a rate of cells budding off from the tail bud. 
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Figure legends 



Figure 1: Top, left: A diagram of the Writ model with a feed-back loop over Dkkl. 
Included are only members of the Wnt pathway that are important for the under- 
standing of the negative feedback loop. When the Wnt signal is on, then Axin gets 
degraded at the LRP5/6 complex and /3-catenin can act as a transcription factor of 
the Wnt inhibitor Dkkl. Vice- versa, when the Wnt signal is off, due to inhibition 
by Dkkl, then /3-catenin gets degraded. Top,right Simulated time series for a 
selection of variables from the model listed in the bottom panel. The total level 
of Axin (magenta) is low, which complies with the findings of ?. Bottom: The 
terms describing the time dependence of the variables. 



Figure 2: Top: Time series for the Dkkl concentration. Space is introduced 
by letting a cell bud off from the tail bud every 10th min. Thus the cells move 
relatively in the PSM. At the determination front the oscillations arrest. Bottom: 
The level of Dkkl at the determination front. (A) R = 10 min. (B) R = 10 min. 
(C) R = 10 min and the cells have the same initial level of Dkkl TB: tail bud. Sf. 
Somite i, where So is the newly formed somite. 



Figure 3: The parameter c^ LGj ^ is changed, and the amplitude (solid) and the 
period (dashed) are calculated. The green vertical line denotes the value of c^qj^ 
in the reference state. 
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Figure 1 
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Tables 



Table I: Parameters in our model of the Wnt system and their default values 



Parameter 


Process 


Default Value 


K c 


Dissociation constant C 


8nM 


CbC 


Breaking of C 


7min _1 


a 


Degradation of /3-catenin 


2.2 min -1 


K[GA] 


Dissociation constant [GA] 


1.5 nM 


Cb[GA] 


Breaking of [GA] 


4min~ 1 


K[LGA] 


Dissociation constant [LCM.] 


InM 


Cb[LGA] 


Breaking of [LGA] 


10 min -1 


V 


Degradation of Axin 


3.8 mm" 1 


K [LD] 


Dissociation constant [LD] 


0.5 nM 


Cb[LD] 


Breaking of [LD] 


0.02 min" 1 


S B 


Constant source of /3-catenin 


1 nM/min 


Sa 


Constant source of Axin 


0.02 nM/min 


Ctsl 


Transcription of dkkl 


0.02 min" 1 


C-tsc 


Translation of Dkkl mRNA 


0.025 (Nm 2 min) -1 


T Dm 


Average lifetime of dkkl mRNA 


8 min 


td 


Average lifetime of Dkkl 


16 min 


GSK3(3 tot 


Total G level 


45 nM 


Ltot 


Total L level 


15 nM 
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